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We study the process e e — > n n n n 7, with a hard photon radiated from the initial state. 
About 60,000 fully reconstructed events have been selected from 89 ftT 1 of BABAR data. The 
invariant mass of the hadronic final state defines the effective e + c~ center-of-mass energy, so that 
these data can be compared with the corresponding direct e + e~ measurements. From the 47r-mass 
spectrum, the cross section for the process e + e~ — > tv + tv~-k + -k~ is measured for center-of-mass 
energies from 0.6 to 4.5 GeV. The uncertainty in the cross section measurement is typically 5%. 
We also measure the cross sections for the final states K + K~tt + -k~ and K + K~ K + K~ . We observe 
the J/tp in all three final states and measure the corresponding branching fractions. We search for 
V(3872) in J/y>(^ )vr + 7r _ and obtain an upper limit on the product of the e + e~ width of the 

X(3872) and the branching fraction for X(3872) -> J/^tt+tt'. 

PACS numbers: 13.66.Bc, 14.40.Cs, 13.25.Gv, 13.25.Jx, 13.20.Jf 



I. INTRODUCTION 

The idea of utilizing initial-state radiation (ISR) from 
a high-mass state to explore electron-positron processes 
at all energies below that state was outlined in Ref. [1]. 
The possibility of exploiting such processes in high lu- 
minosity (f>- and B-factories was discussed in Refs. [2-4] 
and motivates the study described in this paper. This 
is of particular interest because of the small discrepancy 
between the measured muon g — 2 value and that pre- 
dicted by the Standard Model [5], where hadronic loop 
contributions are obtained from e + e~ experiments at low 
center-of-mass (cm.) energies. The study of ISR events 
at B-factories provides independent and contiguous mea- 
surements of hadronic cross sections in this energy region 
and also contributes to the investigation of low-mass res- 
onance spectroscopy. 

The ISR cross section for a particular hadronic final 
state / (excluding the radiated photon) is related to the 
corresponding e+e~ cross section 07 (s) by: 



Events corresponding to 



7 provide the 



d<jf(s, x) 
dx 



W(s,x)-(Tf(s(l-x)) 



(1) 



where x = 2E 7 /y/s; E 1 is the energy of the ISR photon 
in the nominal e + e~ cm. frame; is the nominal e + e~ 
cm. energy; and \J s(\ — x) is the effective cm. energy 
at which the final state / is produced. The function 



W(s,x) = 0- (1 + 6) -x 



1 + 



(2) 



(see for example Ref. [4]) describes the energy spectrum 
of the ISR photons, where — 2a/-K ■ (21n( v / s/m e ) — 1) 
and S takes into account vertex and self-energy correc- 
tions. At the r(4S) energy, 10.58 GeV, = 0.088 and 
S = 0.067. ISR photons are produced at all angles. For 
the present study it is required that the hard ISR photon 
be detected in the electromagnetic calorimeter (EMC) of 
the BABAR detector. Our acceptance for such photons is 
10-15% [4]. 



*Also with Universita dclla Basilicata, Potenza, Italy 
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ISR-luminosity normalization for the hadronic cross sec- 
tion measurements. For a hadronic final state /, the 
Born cross section at center-of-mass energy squared s', 
Uf(s'), is obtained by relating the observed number of 
events, dNf 7 , in an interval ds' centered at s', to the cor- 
responding number of radiative di-muon events, dN^ 7 , 
by means of 



dN^ ■€/■(! + 4sr) 



Ce+c 



^V-(s') • (3) 



Here s' = s(l — x), e w and e/ are detection efficiencies, 

and (1+^fsr) an< ^ are corrections for the possi- 

bility that the detected hard photon may be the result of 
final-state radiation (FSR). This correction is important 
for di-muon events, but is negligible for most hadronic 
final states. The Born cross section er e + e -^+p- ( s ') is 
used. The radiative corrections to the initial state, the 
acceptance for the ISR photon, and the virtual photon 
properties are the same for yU + ^i~ and /, and cancel in 
the ratio. 

An important advantage of ISR data is that the entire 
range of effective cm. energies is scanned in one experi- 
ment. This avoids the relative normalization uncertain- 
ties that inevitably arise when data from different exper- 
iments, or from different machine settings, arc combined. 

A disadvantage of the ISR measurement is that the 
mass resolution is much poorer than can be obtained in 
a direct annihilation. The resolution and absolute energy 
scale can be monitored directly by the width and mass 
of the J/ip resonance produced in the reaction e + e~ — > 
J /tin- By using a kinematic fit to this reaction, we find 
the resolution to be about 8 MeV/c 2 for decays of J/t/j 
in the ji + mode [7]. 

Preliminary studies of e + e~ — > and some mul- 

tihadron ISR processes have been performed with BABAR 
data [6-8]. These demonstrated good detector efficiency 
and particle identification capability for events of this 
kind. 

This paper reports analyses of the 7r + 7r~7r + 7r~, 
K + K ~n + ir~ and K + K~ K + K~ final states produced in 
conjunction with a hard photon, assumed to result from 
ISR. While BABAR data arc available at effective cm. en- 
ergies up to 10.58 GeV, the present analysis is restricted 
to energies below 4.5 GeV because of backgrounds from 
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Y(4S) decays. A clear J/ip signal is observed for each of 
these hadronic states and the corresponding J/ip branch- 
ing fractions are measured. A search for the A (3872)7 
process with the A (3872) decay to J/ipir + Tr~ is also car- 
ried out. 



II. THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II asymmetric e + e~ storage 
ring. The total integrated luminosity used is 89 fb _1 , 
which includes data collected at the T(4S) resonance 
mass (80 fb _1 ), and at cm. energy 40 MeV lower 
(9 fb- 1 ). 

The BABAR detector is described elsewhere [9] . Final 
states with four charged particles are reconstructed in the 
BABAR tracking system, which comprises the Silicon Ver- 
tex Tracker (SVT) and the drift chamber (DCH). Separa- 
tion of pions and kaons is accomplished by means of the 
Detector of Internally Reflected Cherenkov Light (DIRC) 
and energy-loss measurements in the SVT and DCH. The 
hard photon is detected in the EMC. Muon identification 
is provided by the Instrumented Flux Return (IFR) and 
this information is used to select the /i + /i~7 final state 
and in the A (3872) search. 

The initial selection of candidate events requires that 
a high-energy photon in the event with E^ m > 3 GeV 
be found recoiling against four good quality charged 
tracks with zero net charge. Events having such a high- 
energy photon together with an odd number (> 3) of 
good charged tracks are also selected, for the purpose 
of making estimates of the tracking efficiency. Each 
charged track is required to originate close to the in- 
teraction region, to have transverse momentum greater 
than 0.1 GeV/c and to have a polar angle in the labora- 
tory frame with respect to the collision axis in the range 
from 0.4 to 2.45 radians. These selections guarantee the 
quality of the charged tracks in the DCH. Events with 
electrons and positrons are removed on the basis of asso- 
ciated EMC-energy deposition and energy-loss (dE/dx) 
information from the DCH. 

In order to study the detector acceptance and effi- 
ciency, we developed a special package of simulation pro- 
grams for radiative processes. The simulation of the 
7r + 7r~7r + 7r~7 final state is based on the generator de- 
veloped by Kuehn and Czyz [10]. The model assumes 
ai(1260)7r dominance [18, 19], so that many of the events 
contain a pair of pions from a p meson due to the decay 
ai(1260) — > p°ir. No corresponding generator exists for 
the K + K~ir + ir~ and K + K~ K + K~ final states, so these 
reactions were simulated according to phase space. 

Multiple soft-photon emission from the initial-state 
charged particles is implemented with the structure- 
function technique [11, 12], while extra photon radiation 
from the final-state particles is simulated by means of 
the PHOTOS package [13]. The accuracy of the radia- 
tive corrections is about 1%. 



A sample of about 400k events were generated with 
these tools and passed through the detector response 
simulation [14]. These events were then reconstructed 
through the same software chain as the experimental 
data. Variations in detector and background conditions 
were taken into account. 

For purposes of background estimation, a large sam- 
ple of events from the main ISR processes (27r7, 37r7 ... 
67T7, 2A7T7 ...) was simulated. This sample exceeded 
the expected number of events in the dataset by a fac- 
tor of about three. In addition, the expected numbers 
of e + e _ — > qq (q = u,d,s,c) events were generated via 
JETSET [15] and e+er t+t~ via KORALB [16] in or- 
der to estimate non-ISR-type background contributions. 
The cross sections for the above processes are known with 
about 10% accuracy or better, which is sufficient for the 
background contribution study. 



III. THE KINEMATIC FIT PROCEDURE 

The initial sample of candidate events is subjected to 
a constrained kinematic fit in conjunction with charged- 
particle identification to extract events corresponding to 
the final states of interest. 

For each particular four-charged-particle candidate, 
and for each possible combination of particle types (i.e. 
4tt, 2K2ir or 4K), a one-constraint kinematic fit is per- 
formed without using information from the detected pho- 
ton candidate. Because of the excellent resolution of 
the DCH, the three-momentum vector of the photon is 
better determined through momentum conservation than 
through measurement in the EMC. As a consequence, the 
calibration accuracy of the EMC and its alignment with 
respect to the DCH do not contribute to the systematic 
uncertainties. The initial e + e~ and final-state charged- 
particle four-momenta and their covariance matrices are 
taken into account. The momentum vector of the pho- 
ton reconstructed by the fit in the laboratory frame is 
required to have polar angle 9^ in the range from 0.35 
to 2.4 radians and to match the measured polar angle 
^mcas Q £ a canc jidate photon in the EMC within 50 mrad. 

The corresponding azimuthal angles, effi and </>™ eas , are 
also required to agree to this same tolerance. These an- 
gular criteria reduce the background by a factor of about 
two with no noticeable loss of signal. Finally, the polar 
angle 0^ °f each charged track after the fit has to sat- 
isfy 0.45 < 9^ < 2.4 radians in order to fall within the 
acceptance of the DIRC, which provides about 80% kaon 
identification efficiency. 

The fit for the four-pion final-state hypothesis is re- 
tained for every event. If only one track is identified as 
a kaon, or if two oppositely-charged kaons are identified, 
the K + K~tt + it~ fit is also retained. Finally, if two, three 
or four kaons are identified, the four-kaon fit is applied. 

For events with only three charged tracks recoil- 
ing against a candidate photon, the measured four- 
momentum vector of the photon and its covariance ma- 
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80 . 100 
X(4Jt) 

FIG. 1: The one-constraint \ 2 distributions for data (up- 
per histogram) and MC simulation (shaded histogram) four- 
charged-track events fitted to the four-pion hypothesis. The 
cross-hatched histogram is the estimated background contri- 
bution from non-ISR events obtained from JETSET. The sig- 
nal and control regions are indicated. 




3 4 
m(JtVjtV) (GeV/c ) 

FIG. 2: The four-pion invariant mass distribution for the sig- 
nal region of Fig. 1. The points indicate the background esti- 
mated from the difference between data and MC events for the 
control region of Fig. 1, normalized to the difference between 
data and MC events in the signal region of Fig. 1. The cross- 
hatched histogram corresponds to the non-ISR background of 
Fig. 1. 



trix are used in a one-constraint kinematic fit which as- 
sumes, as appropriate, that only a charged pion or kaon 
is undetected. These events are used in the efficiency 
studies described below. 



IV. THE tt+tt-tt+tt- FINAL STATE 
A. Additional selection criteria 

The results of the one-constraint fit to the four 
charged-track candidates are used to make the final se- 
lection of the four-pion sample. We require xln < 30 for 
the four-pion hypothesis, and that any accompanying fit 
to the 2K2tt hypothesis have X2K2n > 10- We estimate 
that these requirements reduce the contamination of the 
4.7T sample by 2K2-K events to about 1% at the cost of 
about 2.4% of the signal events. 

The one-constraint-fit x 2 distribution for the four-pion 
candidates is shown as the upper histogram of Fig. 1, 
while the shaded region is for the corresponding MC- 
simulated pure 47r7 events. The experimental distribu- 
tion has a contribution from the background processes 
but the MC-simulated distribution is also much broader 
than the usual one-constraint \ 2 distribution. This is due 
to multiple-soft-photon emission in the initial state and 
radiation from the final-state charged particles, neither 
of which is included in the constrained fit but which exist 
in data and MC simulation. To illustrate the difference 
of distributions on Fig. 1 , the MC-simulated \ 2 distribu- 
tion is normalized to the data in the region % 2 < 1 where 



contamination of the background events and multiple soft 
ISR and FSR is lowest. 

The cross-hatched histogram in Fig. 1 represents the 
non-ISR background contribution obtained from the 
JETSET simulation of quark-antiquark production and 
hadronization and does not exceed 3%. 

The region 30 < xi-n < 60 is chosen as a control region 
for the estimation of background from other ISR and non- 
ISR multi-hadron reactions. The procedure followed is 
described in the next section. 

The signal region of Fig. 1 contains 67,063 data and 
71,210 MC events, while for the control region the corre- 
sponding numbers are 4,887 and 2,820 respectively. 



B. Background estimation 

MC simulation of the t + t~ final state and ISR produc- 
tion of multi-hadron final states other than 7r + 7r~7r + 7r~ 
shows that such states would yield a background in the 
selected four-pion sample that would exhibit a relatively 
flat contribution to the xL- distribution. We subtract the 
shaded histograms of Fig. 1 from the plain one and the 
resulting histogram is well described by MC simulation 
of background processes. The background contribution 
to any distribution other than \ 2 is estimated as the dif- 
ference between the distributions in the relevant quantity 
for data and MC events from the control region of Fig. 1, 
normalized to the difference between the number of data 
and MC events in the signal region. 
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For example, Fig. 2 shows the four-pion invariant mass 
distribution up to 4.5 GcV/c 2 for the signal region of 
Fig. 1. The points with error bars show the ISR back- 
ground contribution obtained in the manner described 
from the control region of Fig. 1. The cross-hatched 
histograms in Fig. 1 and Fig. 2 represent the non-ISR 
background contribution obtained from the JETSET MC 
simulation. The x 2 distribution for the non-ISR events is 
not flat, and we estimate their relative contribution to the 
signal region from the production cross-section and the 
integrated luminosity. Both backgrounds are small at low 
mass, but the non-ISR background accounts for almost 
half of the observed data at approximately 4 GeV/c 2 . The 
data show a strong peak around 1.5 GcV/c 2 followed by 
a shoulder near 1.9 GeV/c 2 . Narrow signals are apparent 
at the J/tp and the tp(2S) masses, although the latter is 
due to ip(2S) — ► Tr + n~J/ip J/ip — > n + n~ with the muons 
being treated as pions. 

Accounting for uncertainties in cross sections for back- 
ground processes and statistical fluctuations in the num- 
ber of simulated events, we estimate that this procedure 
for background subtraction results in a systematic un- 
certainty of less than 1% in the number of signal events 
in the 1-3 GeV/ c 2 region of four-pion mass, but that it 
increases to 3-5% in the region above 3 GeV/c 2 and to 
roughly 10% in the region below 1 GeV/c 2 . 

By selecting a "background-free" 47r7 sample with only 
four charged tracks and only one photon (about 10% of 
events) we can compare x 2 distributions for data and MC 
events up to x 2 =1000. We estimate that for a xlrr < 
30 selection the net signal size should be increased by 
(3 ± 2)% to allow for a slight shape difference between 
the MC and experimental x 2 distributions. 



C. Tracking efficiency 

We measure the track-finding efficiency with events 
that have three charged-particle tracks and a hard pho- 
ton. These events are subjected to a one-constraint fit 
(as described in Sec. Ill above), which yields the three- 
momentum vector of the missing charged pion in the lab- 
oratory frame assuming this is the only undetected track. 
If the chi-squared of the fit is less than 30 and this vec- 
tor lies within the acceptance of the DCH, the event is 
included in the data sample. 

The four-pion mass distribution obtained in this way 
for three-charged-track events is shown in Fig. 3(a), to- 
gether with the ISR (points with errors) and non-ISR 
(cross-hatched histogram) background estimated as de- 
scribed above. The behavior is similar to that observed 
in Fig. 2. This is exhibited explicitly in Fig. 3(b), where 
the ratio of the three- to four-charged track pion mass 
distributions (after background subtraction) is shown as 
a function of four-pion invariant mass for data (open 
points) . When the MC-simulated data are treated in the 
same way, the solid points of Fig. 3(b) are obtained. The 
same absence of mass dependence is observed for data 




m(% + %'% + %) (GeV/c 2 ) 




m(% + %'% + %) (GeV/c 2 ) 

FIG. 3: (a) The four-pion invariant mass distribution ob- 
tained from fits to three-charged-track events in which an ob- 
served photon is used to constrain the unobserved fourth pion; 
ISR and non-ISR background contributions are indicated by 
the points with error bars and the cross-hatched histogram, 
respectively, (b) The four-pion mass dependence of the ratio 
between the three- and four-charged-track distributions for 
data (open circles) and MC simulation (solid squares). 



and MC events, but the MC simulation yields a smaller 
fraction of three-track events than observed for the data. 
This difference is (3.0 ±0.3 ±2.0)%, where the systematic 
error is estimated from the slight difference in the mass 
dependence seen in Fig. 3(b), from the uncertainties in 
background subtraction and from a slight difference in 
angular dependence. This uncertainty increases to about 
10% for the mass region below 1 GeV/c 2 . The systematic 
difference is used to correct the observed signal size for 
the difference in net track-finding efficiency between that 
obtained from MC simulation and that observed in the 
experiment. 



D. Detection efficiency from simulation 

The selection procedures applied to the data are also 
applied to the MC-simulated event sample. The re- 
sulting four-pion invariant-mass distribution is shown in 
Fig. 4(a) for the signal and control (shaded histogram) re- 
gions. There is good qualitative agreement with the mass 
distribution of Fig. 2, except that no attempt was made 
to simulate the J/ip and ip(2S) signals observed in the 
data. The mass dependence of the detection efficiency 
is obtained by dividing the number of reconstructed MC 
events in each 25 McV/c 2 mass interval by the number 
generated in this same interval. The result is shown in 
Fig. 4(b); the curve is obtained from a polynomial fit to 
the distribution. The efficiency increases from 20% at 
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FIG. 4: (a) The four-pion mass distribution from MC simula- 
tion for the signal and control (shaded) regions of Fig. 1. (b) 
The mass dependence of the net reconstruction and selection 
efficiency obtained from simulation. 
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FIG. 5: (a) The distribution in track-pair opening angle for 
the minimum of the six values possible for each event; (b) the 
distribution in polar angle, and (c) the transverse momen- 
tum distribution for all pions from all events. All quantities 
are in the laboratory frame; the points are for data and the 
histograms are obtained from MC simulation. 



about 0.8 GeV/c 2 to a maximum of 40% near 1 GeV/c 2 , 
and thereafter falls off gradually with increasing mass to 
about 26% at 4.5 GeV/c 2 . This efficiency estimate takes 
into account the geometrical acceptance of the detector 
for the final-state photon and the charged pions, the inef- 
ficiency of the several detector subsystems and event-loss 
due to additional soft-photon emission from the initial 
and final states. 

As mentioned in Sec. II, the model used in the MC 
simulation assumes that the four-pion final state results 
predominantly from the ai(1260)7r quasi-two-body pro- 
duction process [10]. A contribution from /o(1370)p(770) 
is incorporated also. In general, this model describes well 



FIG. 6: The angular distribution of the lowest-momentum 
pion (left) and of the three most energetic pions (right) in 
the four-pion rest frame with respect to the direction of the 
four-pion system in the laboratory frame for the five regions 
of four-pion mass indicated. The points are data, and the 
histograms are MC simulation. 



the distributions in many of the kinematic variables char- 
acterizing the four-pion final state. Some examples are 
shown in Figs. 5 and 6, in which the points with error 
bars represent data while the histograms are obtained 
from MC simulation; for both figures, the J/tp and ^(2S) 
regions have been excluded. Figure 5(a) shows the distri- 
bution in V'min, the minimum charged-pion-pair opening 
angle for each event, while Fig. 5(b) and Fig. 5(c) repre- 
sent the distribution in polar angle, # c h, and transverse 
momentum, pt, respectively, for all final-state pions. All 
quantities are calculated in the laboratory frame, and 
the overall agreement between MC simulation and data 
is very good. Figure 6 compares the distributions in cos 9, 
where 9 is the angle between a charged pion in the four- 
pion rest frame, and the direction of the four-pion system 
in the laboratory frame. The distributions are presented 
for the five regions of four-pion mass indicated (mass in- 
creasing from top to bottom); the left column is for the 
lowest-momentum pion in the four-pion rest frame, and 
the right column sums the distributions for the others. 
Data and MC are in relatively good agreement up to 
about 2 GeV/c 2 , but but above this value some small 
discrepancies appear. 

In the four-pion rest frame, the angular acceptance is 
rather uniform. Changes in the ai(1260) and /o(1370) 
resonance parameter values within the ranges of their 
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uncertainties produce little effect. However, simulation 
without resonances using only four-pion phase space does 
produce discernible deviations from the observed angu- 
lar distributions, and changes the overall acceptance by 
about 2%. This value is taken as an estimate of system- 
atic uncertainty in the acceptance associated with the 
simulation model used. 
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FIG. 7: The integrated ISR luminosity per 0.1 GeV in effec- 
tive cm. energy obtained from e + e + — > h + /j,~~/ events. The 
curve represents the fit used in the luminosity calculations. 
The point at the p(770) mass position was excluded from the 
fit because of pion misidentification feed-through from the 
7r + 7r~7 final state. The points at the J/lp mass position were 
also excluded because of ISR production of J/ip followed by 
decay to /x + [i~ . 



E. Cross section for e + e 



Data from the reaction e + e~ — * /i + /i _ 7 are used 
to convert the invariant-mass distribution for an ISR- 
produced hadronic final state to the energy dependence 
of the corresponding e + e~ cross section. The invariant 
mass of the muon pair rmT defines an effective e + e~ cm. 



collision energy, E c , m .. The differential luminosity, dC, 
associated with the interval dE cm . centered at effective 
collision energy E c . m . is then obtained from 



dC(E c . m .) 



where E r 



dN^(E c 



(1 



dN„ 



T H+li- (E c .m ) ■ (1 + <5 vac ) 

(4) 

is the number of muon 



pairs in the mass interval dm 



/</< 

inv 



dE c 



is the 



acceptance, corrected for muon identification and soft- 
photon emission; (1 + <5pg R ) corrects for hard photon 
emission from final-state muons; cr„+„- (-E c . m .) is the 
e + e~ — > Born cross section at center-of-mass en- 

ergy £"c.m.; and (1 + <5 va c) is the corresponding vacuum 
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FIG. 8: The e + e~ cm. energy dependence of the 7r + 7r"~7r + 7r~ 
cross section measured with ISR data. The events due to 
"4>(2S) — > n + ir~J/ip with J/tjj — > (see Fig. 2) have been 

removed. Only statistical errors are shown. 



polarization correction [20]. For the /i + /i~7 sample ob- 
tained from a BABAR integrated luminosity of 89 fb -1 , 
the dependence of the resulting differential luminosity 
on E C Tn is shown in Fig. 7 in units of [nb _1 /0.1 GeV]. 
From a detailed study of the e + e~ — > /i+^~y detection 
and identification efficiency described in detail in Ref. [7] 
and comparison of the observed invariant-mass spectrum 
with theoretical calculations, we estimate the systematic 
uncertainty associated with luminosity determination to 
be 3%. 

The four-pion e + e~ cross section can then be calcu- 
lated from 



er(7r + 7r 



~7T + 7r 



)(£c.m.) = 



dN ^l( E c.m.) 



dC{E c . m .) 



^(Elm.) ' 

(5) 

Ec.m. with 



where Eg m = 

4tt 



EJ? r 



— ™4tt — 



"hnv 

tern; 



the invariant mass of the four-charged-pion sys- 
dNin-y is the number of selected four-pion events 
after background subtraction in the interval dE c . nl . and 
e^i. c (i? c . m .) is the corresponding detection efficiency ob- 
tained from the MC simulation. The factor £4°" takes 
into account the difference between the \ 2 distributions 
for data and MC events, and the tracking-efficiency dis- 
crepancies discussed in Sec. IV B and Sec. IV C respec- 
tively. 

Since the four-pion cross section calculation involves 
the ratio of the numbers of observed 47r7 and fifi'y events, 
corrections related to multi-soft-photon emission in the 
initial state and to the detection efficiency of the ISR 
photon cancel, these being the same for both reactions. 
Also, since dC has been corrected for vacuum polariza- 
tion and final-state soft-photon emission, the four-pion 
cross section measured in this way includes effects due to 
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vacuum polarization and final-state soft-photon emission. 

We studied the resolution in four-pion mass with MC 
simulation and found that the r.m.s. deviation varied 
from 6.2 MeV/c 2 at mass 1.5 GeV/c 2 to 7.5 MeV/c 2 at 
3 GeV/c 2 . Since the cross section has no sharp peaks 
(J/ip region is discussed below) and the measurements 
are presented in mass intervals of 25 MeV/c 2 , the resolu- 
tion has negligible effect on the measured energy depen- 
dence. 

The energy dependence of the cross section for the re- 



action e + e 



after all corrections is shown 



in Fig. 8. It reaches a peak value of about 30 nb near 
1.5 GeV, with a shoulder at 1.9-2.1 GeV, followed by 
a monotonic decrease toward higher energies perturbed 
only by a small peak at the J/ip mass position. The 
events due to ip{2S) — > ir + tt~ J j 'ip with J/tfj — > ^+ [i~ , 
seen in Fig. 2, have been removed. The luminosity, num- 
ber of events, and corrected cross section for each 25 McV 
interval are presented in Table III. For g — 2 calculations, 
vacuum polarization contributions should be excluded. 
Suitably modified cross section values are presented in 
the last column of Table III. 
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FIG. 10: A detailed view of the energy dependence of the 
e + e~ — > ir + ir~ -k + ir~ cross section near threshold; All exper- 
imental points are indicated as in Fig. 9. Only statistical 
errors are shown. 



We checked the stability of the measured cross sec- 
tion by comparing the results obtained under conditions 
of different DCH voltage, which affects track-hnding effi- 
ciency, and results obtained from data taken at the T(4S) 
mass and below the BB threshold. No significant dis- 
crepancies were observed. 
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FIG. 9: The energy dependence of the e + e~ — > ir + tt~ ir + ir~ 
cross section obtained with BABAR ISR data (black points) in 
comparison with that from direct e + e~ production measure- 
ments. Only statistical errors are shown. 



F. Summary of systematic studies 

The measured four-pion cross section values shown in 
Fig. 8 and summarized in Table III include only statis- 
tical errors. The systematic errors discussed in previous 
sections are summarized in Table I, along with the cor- 
rections that were applied to the measurements. 

The two systematic corrections applied to the mea- 
sured cross sections sum up to +6% with half of this value 
taken as a systematic uncertainty. The systematic errors 
that cancel in the ratio to /X//7 are the photon detection 
efficiency and the ISR soft-photon radiative correction 
uncertainty. 

G. Physics results 

The four-charged-pion cross section measured by 
BABAR can be compared with existing e + e _ measure- 
ments only up to 2.0 GeV — the maximum cm. energy 
up to which measurements of this channel have been 
published. Figures 9 and 10 show the cross section in 
comparison with all existing e + e _ data for cm. energies 
in the 0.7-2.2 GeV range. 

The measured cross section is in good agreement with 
the precision data taken at VEPP-2M by SND [21] and 
CMD-2 [18, 22] in the energy range 0.7-1.4 GeV, as well 
as with data obtained at DCI by DM2 [23] in 1.4-2.0 GeV 
range. The systematic errors for BABAR data are compa- 
rable to, or smaller than, those for these experiments. 

Different mass combinations were studied in data and 
MC events to search for any structures or states not in- 
cluded in the simulation. Figure 11 shows the scatter- 
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TABLE I: Summary of systematic errors for e + e — > tt + tv tt + tt cross section 



Source 


Correction applied 


Systematic error 


Luminosity from /xprf 




3% 

5% for m i7T < 1.0 GeV/c 2 


MC-data difference in x 2 < 30 requirement 


+3% 


2% 


Background subtraction 


- 


1% 


A/lli-H at" a rl 1 ttpvpi 1 pp l ti p k 1 Pit;t;pt; 


+3% 


10% for rruvr < 1.0 GeV/c 2 
3% for m 47r > 3.0 GeV/c 2 
2% 


Radiative corrections accuracy 
Acceptance from MC (model-dependent) 




1% 

2% for < 3 GeV 
15% for the rest 


Total 


+6% 


12% for mi-,, < 1.0 GeV/c 2 


(assuming no correlations) 




5% for 1.0 < m i7I < 3.0 GeV/c 2 
16% for m 4 ^ > 3.0 GeV/c 2 



DATA MC 




m(7t + 7tV"7f ) (GeV/c 2 ) m(7i + 7i _ 7i + 7i ) (GeV/c 2 ) 



FIG. 11: The two-pion (bottom) and three-pion (top) vs. 
four-pion invariant mass distributions for data (left) and sim- 
ulation (right); signals corresponding to J jip and ^p(2S) pro- 
duction are present in the data, but are not included in the 
simulation. Otherwise, agreement between data and MC sim- 
ulation is quite good. 



plots of 37T- and 27r-mass versus 47r-mass for data and 
MC events. Good agreement is seen except for narrow 
regions around the J/tp and ip(2S) masses, which are not 
included in the simulation. 

In order to make a more detailed study, five intervals 
of 47r-mass are selected: (1) 1.0-1.4 GeV/c 2 (for com- 
parison with CMD-2), (2) 1.4-1.8 GeV/c 2 (peak in cross 
section, see Fig. 8), (3) 1.8-2.3 GeV/c 2 (shoulder), (4) 
2.3-3.0 GeV/c 2 , and (5) 3.0-4.5 GeV/c 2 (narrow regions 
around the J ftp and ip(2S) are excluded). Figure 12 
shows the one-dimensional projections of Fig. 11 of the 
27r- and 37r-mass for the above five regions for data and 



MC events. In these distributions we subtract the back- 
ground using control samples in the x 2 distributions from 
data and JETSET simulation, as described above. 

There is excellent agreement between data and MC 
events in the 1.0-1.4 GeV/c 2 region where ai(1260) pro- 
duction is severely limited by the available phase space, 
and /o(1370) production is almost entirely excluded. Dis- 
crepancies begin to appear in the higher mass regions, 
and in particular a relatively narrow bump in the ir + ir~ 
combinations at about 1.3 GeV/c 2 for the 2.3-3.0 GeV/c 2 
region is not reproduced by the simulation. The world 
averages reported by the Particle Data Group (PDG) [17] 
for the ai(1260) mass and width are not well-determined; 
in the simulation we use 1.33 GeV/c 2 and 0.57 GeV re- 
spectively. These values were obtained from a combined 
analysis of CLEO and CMD-2 data [19]. 

The 37r-mass data of Fig. 12 (1.8-2.3 GeV/c 2 ) seem to 
favor a lower ai(1260) mass value than the 1.33 GeV/c 2 
used in simulating the ai(1260). In the simulation, 
the /o(1370) mass and width were set to 1.3 GeV/c 2 
and 0.6 GeV respectively. The 27r-mass data for the 
2.3-3.0 GeV/c 2 region seem to indicate the need for a 
significantly smaller width, although the peak around 
1.3 GeV/c 2 may also be due to significantly higher pro- 
duction of /2(1270) than in the present model. 

In Fig. 12, the right-hand column of plots corresponds 
to the highest di-pion mass value when one other di-pion 
mass is within 25 MeV/c 2 of the p mass. The fiat, broad 
band below 1 GeV/ c 2 is the reflection of the p band in 
ai(1260) decay. The comparatively narrow peak around 
1.3 GeV/c 2 in the 2.3-3.0 GeV/c 2 region indicates the 
need for /o(1370)p in the simulation. 

A full partial wave analysis is required in order to arrive 
at a more precise interpretation of the data. However, 
this requires a simultaneous analysis of the 7r + 7r~7r 7r° 
final state, which is beyond the scope of this study. 
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FIG. 12: The three-pion (left) and 7r + 7r~ (center) mass distributions; the points represent data, the histograms simulation, and 
there are four entries per event. Five regions of four-pion mass are indicated in the center plots. For events having a tt + tt~ 
mass combination in the p(770) region, the distribution of the highest mass other 7r + 7r~ pair is shown in the plots to the right. 



V. THE K+R- 



FINAL STATE 



The constrained fit of the four-charged-track events to 
the hypothesis of two oppositely charged kaons and two 
charged pions, where at least one of the kaons has pos- 
itive particle identification, allows us to select this final 
state. Figure 13 shows the x 2 distributions for both data 
and simulation, where the simulation of the K + K~ir + ir~ 
reaction uses a point-like matrix element with a cross sec- 
tion energy dependence close to that which we observe 
experimentally, and all radiative processes are included. 
Also shown is the estimated contribution from the 47r 
final state arising from particle misidentification, as ob- 
tained from simulation. 

Figure 14(a) presents the simulated mass distribution 
for the 2K2-K events; the mass dependence of the effi- 
ciency, calculated as a ratio of selected to generated 2K2n 
MC events, is shown in Fig. 14(b). 

The selections xIk2-k < 20 > xL > 30 and x\k > 20 
are used, leaving a negligible number of AK and 0.5% 
of 47r events (from simulation) in the final sample. The 
background subtraction procedure which uses the con- 
trol sample with 20 < X2K2n < 40 is applied. The back- 
ground from 47T events is not subtracted by this proce- 



dure, and 0.5% of the events shown in Fig. 2 are used to 
make an additional correction. 

In Fig. 15 we show the K + K~ tt + tt~ invariant mass 
distribution for the events that passed the selection pro- 
cedure. The points with error bars show the background 
distribution from the control sample with 20 < x\k2-k < - 
40 and the shaded histogram is the expected non-ISR 
background from JETSET MC simulation. Both are used 
for background subtraction. A very clear J ftp signal is 
seen. 

Using the number of observed events, efficiency, and 
ISR luminosity, we obtain the e + e~ — > K + K~7t + Tr~ 
cross section shown in Fig. 16, which also displays the 
DM1 data [25]. Table IV presents the cross section in 
25 MeV bins, together with the "undressed" (no vacuum 
polarization) cross section. The systematic errors are 
dominated by uncertainty in the acceptance simulation 
(10%) and the difference between the kaon identification 
efficiencies for data and MC events (up to 5% per track), 
and are estimated to be 15%. 

Figure 17 shows the Kn mass combinations. Produc- 
tion of A"7r pairs is dominated by the A"* (892) clearly 
seen in Fig. 17(a). There also appears to be evidence of 
K% (1430) production in the 2A27T sample seen in the pro- 
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FIG. 13: The one-constraint x 2 distributions for the four- 
charged-track data events and K + K~ tt + tt~ Monte Carlo 
events (shaded histogram) fitted to the K + K~ tt + tt~ hy- 
pothesis. At least one kaon must be identified. The cross- 
hatched histogram is the estimated background contribution 
from four-charged-pion ISR events. The signal and control 
regions are indicated. 
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FIG. 15: The A' + A'~7r + 7r~ invariant mass distribution for 
the signal region of Fig. 13. The points indicate the back- 
ground estimated from the difference between data and MC 
simulation for the control region of Fig. 13 normalized to the 
difference between the number of data and MC events in the 
signal region of Fig. 13. The shaded histogram corresponds 
to the non-ISR background estimated using JETSET. 
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FIG. 14: (a) The K + K~ir + ir~ mass distribution from simula- 
tion for the signal and control (shaded) regions of Fig. 13; (b) 
the mass dependence of the net reconstruction and selection 
efficiency obtained from simulation. 




FIG. 16: The cm. energy dependence of the e e~ — > 
K + K~ n + n~ cross section obtained from ISR events at 
BABAR compared with the only direct e + e~ measurement, 
by DM1. Only statistical errors are shown. 



jection plot of Fig. 17(a) shown in Fig. 17(c). No struc- 
tures in m(isr + 7r + ) or in m(K~ir~) are seen in Fig. 17(b). 
The Kit mass distribution for the other Kn combination 
for events in the K*°(892) bands of Fig. 17(a) is shown in 
Fig. 17(d); events in the overlap region of the two bands 
are included only once to avoid double-counting. The ab- 



sence of a clear if*°(892)_or if 2 * (1430) signal in Fig. 17(d) 
indicates that J ftr* (892) J FsT* (892) and #*° (892)1^(1430) 
quasi-two-body production reactions are small. 

When events in the K*°(892) bands of Fig. 17(a) are 
removed, the scatter-plot m(7r + 7r~) vs. m(K + K~) in 
Fig. 18(a) shows the presence of the p° and <f> resonances. 
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The 7r + 7i- (Fig. 18(b)) and K + K~ (Fig. 18(c)) mass pro- 
jections from Fig. 18(a) exhibit clear p° and 4> signals, 
respectively. Figure 18(d) shows the 7r + 7r _ mass distri- 
bution for events with m(K + K~) in the <j> region. The 
absence of any p signal is consistent with C conservation; 
there might be slight evidence of /o(980) production. 

The three-body mass combinations are also of poten- 
tial interest. We consider these in two categories. For 
the first category, we require that there be a K^tt^ com- 
bination in one of the K*°(892) bands of Fig. 17(a). The 
K*°(892)K ± invariant mass behavior for such events is 
shown in Figs. 19(a,b). In Fig. 19(b), there is a strong 
peak just below 1.5 GcV/c 2 followed by a shoulder around 
1.6-1.7 GeV/c 2 and a rapid drop in the 1.8-2.0 GeV/c 2 
region. The peak could correspond to a presently un- 
known isovector state decaying to K*°K + K*°K. The 
corresponding K*°ir behavior is shown in Figs. 19(c,d). 
The mass projection in Fig. 19(c) shows a broad struc- 
ture in the region of the ifi(1270) and ifi(1410), both of 
which arc known to decay through K*°ir. The scatter- 
plot ofm(K*°Tr) vs. m(K*°K) in Fig. 19(d) shows that 
the low-mass enhancements of Figs. 19(b,c) are highly- 
correlated. 

For the second category, we exclude events with at least 
one Kit combination in a if* (892) band of Fig. 17(a). 
As shown in Fig. 18(c), a subset of the remaining events 
contains a (j> signal. When the events in the <f> region 
are combined with the remaining ^ , the (fyir^ mass dis- 
tribution of Fig. 20(a) is obtained. This shows a peak 
in the mass region 1.25-1.4 GeV/c 2 , followed by a sec- 
ond broad bump in the 1.5-1.8 GeV/c 2 region. If the 
4> region is excluded, the rather featureless mass distri- 
bution of Fig. 20(b) is obtained. However, the scatter- 
plot of m(7r + 7r - ) vs. m(_fiT : ' = 7r + 7r _ ) in Fig. 20(c) for these 
events shows clear evidence for p production correlated 
with K ± TT + ir~ mass in the 1.2-1.5 GeV/c 2 region. The 
i^i (1270) couples quite strongly to Kp [17], and so may 
be produced in this final state. However, the K i (1410) 
is almost decoupled from Kp [17], and so cannot be the 
source of the events in the 1.4-1.5 GeV/c 2 region. 

There is evidence in the K + K~-k + -k~ final state for 
several interesting, but complex, structures. Detailed un- 
derstanding will require data on K + K~Tr a ir° and also on 
four-body final states involving neutral kaons. 



VI. THE K+K'K+K- FINAL STATE 

The one-constraint fit for the four-charged-kaon hy- 
pothesis with xIk < 20 selection gives a very pure sam- 
ple of this final state. The main background contribu- 
tion is from the reaction with two charged kaons and 
two misidentified charged pions. We reduced this back- 
ground greatly by adding to the initial selection criterion 
the requirements that three or four of the charged parti- 
cles have good kaon identification and that xIk2tt > 20 
for the fit to the reaction with two kaons and two pions. 
Background from the 4ir final state is negligible. 
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FIG. 17: For the K + K~n + ir~ data sample: (a) the scatter- 
plot of the K + ix~ and K~ix + invariant mass values; (b) 
the scatter-plot of the K + ir + and K~tt~ invariant mass val- 
ues; (c) the K + ir~ or K~ir + mass projection of (a); (d) the 
mass distribution for other Kir combinations for events in the 
A"* (892) bands of (a). 
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FIG. 18: For the K + K~-k + -k~ data sample, after removing 
events contributing to the A"* (892) mass regions indicated 
in Fig. 17(a): (a) the scatter-plot of -k + iv~ and K + K~ mass 
values; (b) the -k + it~ mass projection of (a); (c) the K + K~ 
mass projection of (a); (d) the -k + -k~ mass projection of (a) 
for events in the 0(1020) mass peak. 



Figure 21 shows the x 2 distributions for data and simu- 
lation. The simulation of the AK final state uses a point- 
like matrix element with cross section behavior close to 
what we observe experimentally, and with all radiative 
processes included. Also shown is the contribution from 
the 2K2w final state obtained from the simulation (cross- 
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FIG. 19: The 2K2ty events from the K*°(892) bands in 
Fig. 17(a): (a) the K*°(892)K ± vs. K+R--K+-R- invari- 
ant mass distribution scatter-plot; (b) the K*° (892) K mass 
projection; (c) the _K'* (892)7r ± mass distribution; (d) the 
m{K*°TY T ) vs. m(K*°K±) scatter-plot. 
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FIG. 20: The 2K2n events with the if* (892) bands in 
Fig. 17(a) excluded: (a) the K + K~ mass distribution with 
the K + K~ mass in the <j> region; (b) the K~ "7r mass dis- 
tribution with K + K~ in the cj> mass region excluded; (c) the 
7r + 7r~ vs. K ± tt + it~ scatter-plot, and (d) the K ± tt + tt~ pro- 
jection, for the events from (b). 
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FIG. 21: The one-constraint x 2 distributions for four-chargcd- 
track data events and K + K~ K + K~ Monte Carlo events 
(shaded histogram) fitted to the K + K~ K + K~ hypothesis. 
At least three kaons must be identified. The cross-hatched 
histogram is the estimated background contribution from 
K + K~ tt + tt~ ISR events. The signal and control regions are 
indicated. 




m(K + K"K + K") (GeV/c ) 



> 
a 



0.3 



0.2 



0.1 




2 3 4 . 

m(K K~K K~) (GeV/c 2 ) 

FIG. 22: (a) The K + K~K + K " mass distribution from sim- 
ulation for the signal and control (shaded) regions of Fig. 21; 
(b) the mass dependence of the net reconstruction and selec- 
tion efficiency obtained from simulation. 



hatched). 

Figure 22(a) shows the simulated mass distribution for 
AK events, and Fig. 22(b) presents the detection effi- 
ciency calculated as the ratio of selected to generated 
AK MC events. 

In Fig. 23 the K + K~K + K~ invariant mass distribu- 
tion is shown for the data in the signal region, and also for 



background estimated from the 20 < x\k < ^0 control 
sample (points). The shaded histogram shows the ex- 
pected non-ISR background from the JETSET MC sim- 
ulation. A clear signal due to the J/ip is seen. The 
number of signal events as a function of four-kaon mass 
is obtained by subtracting the background contributions 
estimated from the control sample and from the JETSET 
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FIG. 23: The K+ K~ K+ K~ invariant-mass distribution for 
the signal region of Fig. 21. The points indicate the back- 
ground estimated from the difference between data and MC 
simulation for the control region of Fig. 21, normalized to the 
difference between the number of data and MC events in the 
signal region of Fig. 21. The shaded histogram corresponds 
to the non-ISR background estimated from JETSET. 
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FIG. 24: The cm. energy dependence of the e + e~ — > 
K + K~ K + K~ cross section obtained from ISR events at 
BABAR. Only statistical errors are shown. 



MC simulation. 

Using the number of four-charged-kaon events ob- 
served, the acceptance discussed previously, and the ISR 
luminosity, we calculate the cross section for the reaction 
e + e~ — > K + K~K + K~ . The results with statistical er- 
rors only are displayed in Fig. 24, and the values are listed 
in Table V. There are no published electron-positron 
data for comparison. Systematic errors are dominant 
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FIG. 25: (a) The K + K~ and (b) the K ± K + K~ invariant 
mass distributions for events from the K + K~ K + K~ sample. 



for these measurements. These uncertainties are due to 
the absence of a detailed model for the acceptance sim- 
ulation, to uncertainties in the background subtractions, 
and to differences in kaon identification efficiency in the 
MC simulation and the data. The overall estimated sys- 
tematic error is about 25%. The K+R- and K ± K + K~ 
invariant mass distributions are shown in Fig. 25 (four 
entries per event). If there is any <p production it is 
small, and neither mass distribution shows evidence of 
significant structure. 



VII. THE J/tp REGION 



Figure 26 shows an expanded view of the J/tp mass re- 
gion in Fig. 2 for the four-pion data sample. The signals 
from J/tp — > 7T + 7T~7r + 7r~ and tp(2S) — > J/tpir + n~ — > 
^ + /i~7r + 7r~ (without muon identification) are clearly 
seen. The observation of J/tp decaying into four pi- 
ons gives a direct measurement of the 47r-mass resolution 
and the absolute energy scale. A fit with a Gaussian for 
the J/tp peak and a polynomial function for the contin- 
uum, gives (T m (4,j = 8 McV/c 2 and less than 1 MeV/c 2 
difference from the PDG [17] value for the J/tp mass. 
The observed mass resolution agrees with the simulation 
within 10%. 



The observed 270 ± 20 events at 
be used to calculate the branching 
7T + 7r~7r + 7r~. The simulation shows 
diative effects only 89% of the sign 
the Gaussian curve. This value is 
with that obtained from ip(2S) peak 
low) which is used for calculations. 



the J/tp peak can 
fraction for J/tp — > 
that because of ra- 
al events are under 
in good agreement 
87.7±1.3% (see be- 
Using the corrected 
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FIG. 26: The n + ir~ 7v + ir~ mass distribution for ISR- produced 
e + e~ — > 7r + 7r _ 7r + 7r~ events in the J/ip-^)(2S) region; there 
are clear signals at the J/ip and tp(2S) mass positions. 
The shaded region at the latter corresponds to tp(2S) — > 
J/ipir + ir~ , with J/i/> — > fi + n~ and the muons are treated 
as pions. 



number, we can calculate the products: 



N(J/lp — > 7T+7T 7T+7T ) 

dC/dE ■ cmc 
46.8 ±3.5 ±3.3 nbMeV , 
N(J/ip — > 7r + 7r~7r+7r~) • 

6tt 2 • dC/dE-euc • C 
= (1.95 ± 0.14 ± 0.13) x 10 -2 keV 



where 



12790 ± 850 nb- McV 



is the integral over the J/ip excitation curve; T'U^ is the 
electronic width; dC/dE = 25.3 nb _1 /MeV is the ISR 
luminosity at the J/ip mass; cmc = 0.26 ± 0.01 is the 
detection efficiency from simulation with the corrections 
discussed in Sec. IV F; and C = 3.894 x 10 11 nbMeV 2 
is a conversion constant. The systematic error includes 
3% uncertainty from the ISR luminosity, 3% from effi- 
ciency, and 5% uncertainty from background subtraction. 
The subscript "47r" for branching fractions refers to the 
7r + 7T~7r + 7r~ final state exclusively. 

Using ril^ = 5.40±0.17 keV [17], we obtain the result 
Bj/^a* = (3.61 ± 0.26 ± 0.26) x 10~ 3 , substantially 
more precise that the current PDG value Bj^^ 47T = 
(4.0 ± 1.0) x 10~ 3 [17]. The systematic error includes 
3% uncertainty in T'ii^ . 
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FIG. 27: The K+R--K+-W- (top) and K + R-K + K~ (bot- 
tom) mass distributions for ISR-produced events in the J/ip 
region. A clear signal is observed at the J/ip mass. 



The ip(2S) peak corresponds to the decay chain 
ip(2S) — ► J/-07r + 7r~ — » /i + /i~7r + 7r~ with the muons 
treated as pions. The number of events extracted from a 
fit to a Gaussian distribution for the ip{2S) peak and 
a polynomial function for the continuum is 544 ± 27. 
With a radiative correction of 0.89 taken from simula- 
tion, N(ip(2S) -> J/ipTT+n-) = 611 ± 30 events is the 
observed signal. 

The number of ip(2S) events can be obtained with 
much less background if the invariant mass of one 
pair of charged tracks assumed to be muons is within 
±50 MeV/c 2 of the J/ip mass. Events satisfying this 
criterion are shown by the shaded histogram in Fig. 26; 
this ip{2S) peak has 620 ± 25 events after subtraction 
of 10 background events. With this number and the 
544±27 events obtained from the Gaussian fit, the radia- 
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tive correction value can be obtained directly; the value 
is 0.877 ± 0.013, in good agreement with the value 0.89 
from simulation. 

For the ip(2S) we then obtain: 



D 



■ B 



- ■ o-. 



jV(V»(2S) -> J/jjTT+TT-) 

dC/dE ■ e M c 
76.3 ±3.1 ±3.8 nbMeV 



which leads to 



B 



' £>J/d-m.+ ii- ■ 1 I 



N(ip(2S) -» J/iPir+ir 



m 



V>(2S) 



6 • 7T 2 • AC/dE • £mc • C 1 
= (4.50 ± 0.18 ± 0.22) x 10~ 2 keV , 

where dC/dE = 32.4 nb" 1 / McV is the ISR luminosity at 
the ip{2S) mass, and cmc = 0.25 ± 0.01 is the detection 
efficiency from simulation with corrections discussed in 
Sec. IV F. The same systematic errors as for the J/ip 
have been included. 

Using the values rt^ S) = 2.12 ± 0.12 keV and 
Sj/V-m+m- = 0.0588 ± 0.0010 from Ref. [17], we obtain 

S^swam+tt- = 0.361 ± 0.015 ± 0.028 , 

which should be compared to the current world-average 
value, Bv(2S)-J/V"r+*- = 0.317± 0.011 [17]. Almost half 
of our systematic error in B^ 2 s)->J comes from 

the uncertainty in rte 2 ^ ■ 

Alternatively, using the PDG values for the two 
branching fractions, we can extract the e + e~ width of 
the ip(2S), 

r VK2S) = 2.4i ± o.io ± 0.15 keV , 

where the systematic error has been increased due to the 
uncertainty in B^ 2 s)->J/ipir+ir- ■ 

Figure 27 shows the J/ip signal in the K + K~ir + ir~ 
and K + K~ K + K~ modes. The numbers of events under 
the Gaussian curves are 233 ± 19 and 38.5 ± 6.7 respec- 
tively. The mass resolution is about 7 MeV/c 2 for 2K2tt 
and 5 MeV/c 2 for the AK final state. Using the radia- 
tive correction factor 0.877, and €2K2tt — 0.13 ± 0.01 and 
C4K = 0.18 ± 0.02 from simulation, we obtain 

. = (3.36 ± 0.27 ± 0.27) x 10" 2 keV 



B 



Bj/^ak ■ T J J^ = (4.0 ± 0.7 ± 0.6) x 10~ 3 keV . 

The systematic errors are mainly due to the uncertainties 
in acceptance and ISR luminosity. 

Using the PDG value for Y, 
branching fractions 



we calculate the 



B.,/^2K2, = (6.2 ± 0.5 ± 0.5) x 10- 3 
Bj/^k = (7.4 ±1.2 ±1.2) x 10" 4 



to be compared with the current PDG values of (7.2 ± 
2.3)-10~ 3 and (9.2 ± 3.3)-10~ 4 respectively The uncer- 
tainty in Tii^ has been added in quadrature to the sys- 
tematic error estimate. 

Since we have measured the products of o'H^ and 
branching fraction of J/ip decay to Air, 2K2ir, and AK 
it is interesting to compare them with the non-resonant 
cross sections (continuum) at that energy. Using a linear 
approximation of the cross sections from Tables III, IV, 
V around the J/ip peak within ±0.1 GeV (events from 
the peak are excluded), the following cross sections arc 
obtained at the J/ip mass: 

0-4^ = 0.55 ±0.03 nb 
(J2K2-K = 0.48 ± 0.04 nb 
a iK = 0.055 ± 0.009 nb . 

Table II presents the ratios Bji^^f ■ /c e +e-^f 
for / = 47r, 2K2tt, AK. In these ratios all experimen- 
tal systematic errors cancel. Also shown is the ratio 



B, 



j/^h+ij,- -<rUt '/<7e+e-^+»- taken from Ref. [7]. 



TABLE II: Ratios of the J/ip partial production rates to con- 
tinuum cross sections. The result for fi n~ is from Ref. [7]. 



Final state, / 


Bj/^f<r^/a e+e -^ f (MeV) 


7r + 7T _ 7r + 7T _ 


85.1 ±7.9 


K+K-n+TY- 


166 ±19 


K+K-K+R- 


138 ± 32 




84.12 ±0.67 



The ratio obtained for the 4-7T final state is in good 
agreement with that for Indeed, the strong decay 

of the J/ip to An is forbidden by G-parity conservation 
and therefore this decay is expected to be dominated by 
single-photon. No such suppression due to G-parity for 
the strong decay of the J/ip for the other two modes. 
We interpret the significantly larger values of the ratios 
as an indication that the single-photon mechanism is not 
dominant for the J/ip decays to the 2K2n and AK final 
states. 



VIII. THE SEARCH FOR THE X(3872) 

Events in the J/ip mass region can be used to search 
for the new narrow A (3872) state reported by Belle [26]. 
If this state has J PC = 1~~ it could be seen in ISR 
production via the decay to J/ip-K + -K~~ , just as for the 
ip(2S) (see Fig. 26). We acknowledge that such a J PC 
assignment is improbable since it would permit the un- 
observed decay to DD, which would make the A(3872) 
much broader. If the J/ip is selected in its decay to 
two muons, this state can be searched for as a peak at 
3872 MeV/c 2 in the four-charged-track sample. Using 
our ISR data, we can set an upper limit on the product 



B 



X^.J/%l)TT 



+ w - ■ Tf e . To improve background rejection, we 
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FIG. 28: The m(2p2Tv) — m(2/i) difference for events with 
m(2/i) near the J/?/> mass and both particles identified as 
muons. The lines show the ip{2S) and X(3872) selection re- 
gions. 



require positive muon identification on two tracks, in ad- 
dition to having the invariant mass of the two muons fall 
in the interval [3.05-3.15] GeV/c 2 . Figure 28 displays the 
difference between the mass of the four-track final state 
and that of the two muons. The vertical lines show the 
selected regions for those events from tp(2S) decay and 
for the A (3872) search. The number of events within the 
ip(2S) interval is used to estimate the di-muon identifi- 
cation efficiency, which is found to be 0.61 for this study. 
Events in the interval [0.65-1.0] GeV/c 2 are used to es- 
timate the background within the expected region for 
A(3872) decay, namely [0.74-0.82] GeV/c 2 . One event 
is found within this interval where 1.4 are expected from 
background; this yields an upper limit of 3.0 events at 
the 90% confidence level. From this limit an upper limit 
on the product of the branching fraction of the A (3872) 
to J/ipTT + ir~ and the e + e~ width of A(3872) is obtained 
through the chain 



B 



X^>J/4>w+n- 



■ B 



T x = 



N(X -> J/V'tt+tt") -m\ 
6 • 7T 2 • dC/dE ■ e M c ■ C ■ 0.61 



< 0.37 eV , (6) 



where dC/dE = 34.7 nb _1 / MeV is the ISR luminosity 
at the A (3872) mass, and emc = 0.25 ± 0.01 is the ac- 
ceptance from simulation. 

Using the value B J/4 ,^ 2fi = 0.0588 ± 0.0010 [17], we 
extract the upper limit 



D 



X(3872)^7/V7r+7r- 



• < 6.2 eV at 90% C.L. 



This result is the best upper limit to date and can be 
compared with the result of a similar study performed by 
BES [27]. 



IX. SUMMARY 



The photon energy and charged particle momentum res- 
olutions together with the particle identification capa- 
bilities of the BABAR detector permit the reconstruction 
of 7r + 7r~7r + 7r~, K + K~ n + ir~ and K + K~ K + K~ final 
states produced at low effective cm. energy via ISR in 
data taken in the T(4S) mass region. 

The analysis shows that luminosity and efficiency can 
be understood with 2-4% accuracy, and that ISR pro- 
duction yields useful measurements of R, the ratio of 
the hadronic to di-muon cross section values, in the low- 



energy regime of e + e collisions. 



The selected multi-hadronic final states in the broad 
range of accessible energy provide new information 
on hadron spectroscopy. 



The observed e + c 



7r + 7T~7r + 7r~ cross section provides evidence of resonant 
structure, with preferred quasi-two-body production of 
al(1260)7r. For the first time there is an indication of 
a /o(1370)p(770) contribution to the final state. How- 
ever a detailed understanding of the four-pion final 
state requires additional information from states such as 
7r + 7r _ 7r°7r°. 

The cross section measurements for the reaction 
e + e~ — > K + K~tt + tt~ present a significant improvement 
upon existing data with about 15% systematic uncer- 
tainty. In addition, the final state exhibits complex res- 
onance sub-structure. Clear signals for if* (892), 0, 
and p are observed; there is evidence of A,* (1430) pro- 
duction; and there are low-mass enhancements in the 
K*°(892)K ± , A* (892)tt, 0tt, and Kp subsystems. It 
is difficult to disentangle these contributions to the final 
state, and we make no attempt to do so in this paper. 

The energy dependence of the cross section for the re- 



action e + e 



K+K~K+K~ from threshold to 4.5 GeV 



has been measured for the first time with about 20% sys- 
tematic uncertainty. In contrast to the A+ A~7r+7r~ final 
state, the four-kaon state shows no clear mass structure 
in the two- and three-body sub-systems. The absence of 
a clear signal in the K + K~ mass distribution is unex- 
pected. 

The analyzed 89 fbT 1 of BABAR data in the 1.4- 
4.5 GeV/c 2 mass range are already better in quality and 
precision than the direct measurements from the DCI 
and ADONE machines, and do not suffer from the rela- 
tive normalization uncertainties which seem to exist for 
direct measurements of these final states. 

The ISR events allow a study of J/ip and ip(2S) pro- 
duction, and the measurement of the product of decay 
branching fractions and e + e~ width of the J/ip with the 
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best accuracy to date. The results are as follows: 

= (1.95 ± 0.14 ± 0.13) x keV , 

B^,(2S)^j/4>7t+7t- ■ Bj/i[,^2n ■ rfj 2S,) 

= (4.50 ± 0.18 ± 0.22) x 10" 2 keV , 

= (3.36 ± 0.27 ± 0.27) x 10" 2 keV , 

Bj/^^ak ■ 

= (4.0 ±0.7 ±0.6) x 10~ 3 keV . 

The dominance of the singlc-photon-decay mecha- 
nism for J/ip — > 7r + 7r _ 7r + 7r _ has been demonstrated by 
comparison with the continuum cross section e+e~ — > 
7r + 7r~7r + 7r~. 

Under the assumption that the X(3872) has J PC = 
we have obtained the best upper limit to date on 
the product of its branching fraction to J/ipn + ir~ and 
its e + e~ width: 

B X (3S72)^j/^+k- ■ rf e < 6.2 eV at 90% CL. 
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TABLE III: Summary of e + e — ► tv + tt tv + tt cross section measurement. "Dressed" and "Undressed" (without vacuum 
polarization) cross sections are presented. Errors are statistical only. 
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1 
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TABLE IV: Summary of e + e — ► K + K tx + -k cross section measurement. "Dressed" and "Undressed" (without vacuum 
polarization) cross sections are presented. Errors are statistical only. 



Ec. m . (GeV) 


a 


(nb) 


OnoVP 


(nb) 


Ec.m. (GeV) 


a 


(nb) 


OnoVP 


(nb) 


E c . m . (GeV) 


a 


(nb) 


fnoVP 


(nb) 


1.5125 


0.00 


± 


0.00 


0.00 


± 


0.00 


2.5125 


1.27 


± 


0.19 


1.22 


± 


0.18 


3.5125 


0.33 


± 


0.07 


0.32 


± 


0.07 


1.5375 


0.03 


± 


0.05 


0.03 


± 


0.05 


2.5375 


1.20 


± 


0.18 


1.15 


± 


0.17 


3.5375 


0.29 


± 


0.07 


0.28 


± 


0.07 


1.5625 


0.11 


± 


0.05 


0.11 


± 


0.05 


2.5625 


1.58 


± 


0.19 


1.52 


± 


0.18 


3.5625 


0.22 


± 


0.06 


0.21 


± 


0.06 


1.5875 


0.29 


± 


0.08 


0.28 


± 


0.08 


2.5875 


1.22 


± 


0.17 


1.17 


± 


0.16 


3.5875 


0.07 


± 


0.06 


0.07 


± 


0.06 


1.6125 


0.58 


± 


0.13 


0.56 


± 


0.13 


2.6125 


1.02 


± 


0.16 


0.98 


± 


0.15 


3.6125 


0.08 


± 


0.05 


0.08 


± 


0.05 


1.6375 


0.96 


± 


0.15 


0.93 


± 


0.14 


2.6375 


1.03 


± 


0.16 


0.99 


± 


0.15 


3.6375 


0.08 


± 


0.05 


0.08 


± 


0.05 


1.6625 


1.22 


± 


0.19 


1.18 


± 


0.18 


2.6625 


0.83 


± 


0.17 


0.80 


± 


0.16 


3.6625 


0.16 


± 


0.06 


0.16 


± 


0.06 


1.6875 


1.60 


± 


0.22 


1.54 


± 


0.21 


2.6875 


0.69 


± 


0.14 


0.66 


± 


0.13 


3.6875 


0.21 


± 


0.07 


0.17 


± 


0.06 


1.7125 


1.78 


± 


0.25 


1.71 


± 


0.24 


2.7125 


1.18 


± 


0.17 


1.13 


± 


0.16 


3.7125 


0.24 


± 


0.07 


0.22 


± 


0.06 


1.7375 


2.67 


± 


0.30 


2.57 


± 


0.29 


2.7375 


0.67 


± 


0.15 


0.64 


± 


0.14 


3.7375 


0.11 


± 


0.06 


0.10 


± 


0.06 


1.7625 


2.81 


± 


0.33 


2.70 


± 


0.32 


2.7625 


0.98 


± 


0.14 


0.94 


± 


0.13 


3.7625 


0.23 


± 


0.06 


0.22 


± 


0.06 


1.7875 


2.99 


± 


0.34 


2.88 


± 


0.33 


2.7875 


0.78 


± 


0.14 


0.75 


± 


0.13 


3.7875 


0.13 


± 


0.05 


0.12 


± 


0.05 


1.8125 


3.72 


± 


0.36 


3.58 


± 


0.35 


2.8125 


0.72 


± 


0.13 


0.69 


± 


0.13 


3.8125 


0.06 


± 


0.05 


0.06 


± 


0.05 


1.8375 


4.16 


± 


0.40 


4.00 


± 


0.38 


2.8375 


0.77 


± 


0.13 


0.74 


± 


0.13 


3.8375 


0.13 


± 


0.06 


0.12 


± 


0.06 


1.8625 


4.58 


± 


0.41 


4.41 


± 


0.39 


2.8625 


0.63 


± 


0.13 


0.61 


± 


0.13 


3.8625 


0.08 


± 


0.06 


0.08 


± 


0.06 


1.8875 


4.34 


± 


0.41 


4.18 


± 


0.39 


2.8875 


0.28 


± 


0.12 


0.27 


± 


0.12 


3.8875 


0.08 


± 


0.05 


0.08 


± 


0.05 


1.9125 


4.18 


± 


0.38 


4.03 


± 


0.37 


2.9125 


0.48 


± 


0.13 


0.46 


± 


0.13 


3.9125 


0.07 


± 


0.04 


0.07 


± 


0.04 


1.9375 


4.70 


± 


0.42 


4.53 


± 


0.40 


2.9375 


0.53 


± 


0.11 


0.51 


± 


0.11 


3.9375 


0.07 


± 


0.05 


0.07 


± 


0.05 


1.9625 


4.01 


± 


0.38 


3.86 


± 


0.37 


2.9625 


0.44 


± 


0.12 


0.43 


± 


0.12 


3.9625 


0.10 


± 


0.05 


0.10 


± 


0.05 


1.9875 


4.11 


± 


0.38 


3.96 


± 


0.37 


2.9875 


0.54 


± 


0.11 


0.53 


± 


0.11 


3.9875 


0.03 


± 


0.05 


0.03 


± 


0.05 


2.0125 


3.08 


± 


0.33 


2.97 


± 


0.32 


3.0125 


0.48 


± 


0.12 


0.47 


± 


0.12 


4.0125 


0.04 


± 


0.04 


0.04 


± 


0.04 


2.0375 


3.19 


± 


0.33 


3.07 


± 


0.32 


3.0375 


0.39 


± 


0.11 


0.39 


± 


0.11 


4.0375 


0.13 


± 


0.05 


0.12 


± 


0.05 


2.0625 


3.32 


± 


0.35 


3.20 


± 


0.34 


3.0625 


0.46 


± 


0.11 


0.47 


± 


0.11 


4.0625 


0.00 


± 


0.00 


0.00 


± 


0.00 


2.0875 


2.79 


± 


0.32 


2.69 


± 


0.31 


3.0875 


2.40 


± 


0.24 


0.67 


± 


0.07 


4.0875 


0.02 


± 


0.04 


0.02 


± 


0.04 


2.1125 


3.31 


± 


0.32 


3.19 


± 


0.31 


3.1125 


1.30 


± 


0.17 


1.00 


± 


0.13 


4.1125 


0.05 


± 


0.04 


0.05 


± 


0.04 


2.1375 


2.84 


± 


0.32 


2.73 


± 


0.31 


3.1375 


0.59 


± 


0.11 


0.53 


± 


0.10 


4.1375 


0.00 


± 


0.04 


0.00 


± 


0.04 


2.1625 


2.48 


± 


0.28 


2.39 


± 


0.27 


3.1625 


0.48 


± 


0.11 


0.44 


± 


0.10 


4.1625 


0.04 


± 


0.04 


0.04 


± 


0.04 


2.1875 


2.83 


± 


0.30 


2.72 


± 


0.29 


3.1875 


0.32 


± 


0.10 


0.30 


± 


0.09 


4.1875 


0.04 


± 


0.05 


0.04 


± 


0.05 


2.2125 


2.00 


± 


0.25 


1.92 


± 


0.24 


3.2125 


0.57 


± 


0.11 


0.54 


± 


0.10 


4.2125 


0.06 


± 


0.04 


0.06 


± 


0.04 


2.2375 


1.79 


± 


0.24 


1.72 


± 


0.23 


3.2375 


0.46 


± 


0.10 


0.43 


± 


0.09 


4.2375 


0.01 


± 


0.04 


0.01 


± 


0.04 


2.2625 


2.18 


± 


0.26 


2.10 


± 


0.25 


3.2625 


0.40 


± 


0.09 


0.38 


± 


0.08 


4.2625 


0.13 


± 


0.04 


0.12 


± 


0.04 


2.2875 


1.53 


± 


0.23 


1.47 


± 


0.22 


3.2875 


0.32 


± 


0.08 


0.30 


± 


0.08 


4.2875 


0.04 


± 


0.04 


0.04 


± 


0.04 


2.3125 


1.91 


± 


0.24 


1.84 


± 


0.23 


3.3125 


0.40 


± 


0.09 


0.38 


± 


0.09 


4.3125 


0.02 


± 


0.04 


0.02 


± 


0.04 


2.3375 


1.73 


± 


0.24 


1.67 


± 


0.23 


3.3375 


0.25 


± 


0.07 


0.24 


± 


0.07 


4.3375 


0.06 


± 


0.03 


0.06 


± 


0.03 


2.3625 


1.77 


± 


0.21 


1.70 


± 


0.20 


3.3625 


0.40 


± 


0.09 


0.38 


± 


0.09 


4.3625 


0.00 


± 


0.04 


0.00 


± 


0.04 


2.3875 


1.54 


± 


0.22 


1.48 


± 


0.21 


3.3875 


0.24 


± 


0.07 


0.23 


± 


0.07 


4.3875 


0.02 


± 


0.04 


0.02 


± 


0.04 


2.4125 


2.01 


± 


0.23 


1.93 


± 


0.22 


3.4125 


0.22 


± 


0.08 


0.21 


± 


0.08 


4.4125 


0.05 


± 


0.04 


0.05 


± 


0.04 


2.4375 


1.38 


± 


0.20 


1.33 


± 


0.19 


3.4375 


0.00 


± 


0.00 


0.00 


± 


0.00 


4.4375 


0.05 


± 


0.03 


0.05 


± 


0.03 


2.4625 


1.61 


± 


0.22 


1.55 


± 


0.21 


3.4625 


0.24 


± 


0.08 


0.23 


± 


0.08 


4.4625 


0.05 


± 


0.03 


0.05 


± 


0.03 


2.4875 


1.84 




0.21 


1.77 




0.20 


3.4875 


0.21 




0.07 


0.20 




0.07 


4.4875 


0.06 


± 


0.04 


0.06 




0.04 
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TABLE V: Summary of e + e — ► K + K K + K cross section measurement. "Dressed" and "Undressed" (without vacuum 
polarization) cross sections are presented. Errors are statistical only. 



Ec. m . (GeV) 


a (nb) 


fnoVP (nb) 


£cm. (GeV) 


a (nb) 


o"novp (nb) 


Sc.m. (GeV) 


a (nb) 


o"novp (nb) 


2.0312 


0.00 ± 0.01 


0.00 ± 0.00 


2.9062 


0.07 ± 0.02 


0.06 ± 0.02 


3.7188 


0.00 ± 0.01 


0.00 ± 0.01 


2.0938 


0.02 ± 0.01 


0.02 ± 0.01 


2.9688 


0.08 ± 0.02 


0.08 ± 0.02 


3.7812 


0.05 ± 0.01 


0.04 ± 0.01 


2.1562 


0.05 ± 0.01 


0.05 ± 0.01 


3.0312 


0.06 ± 0.02 


0.06 ± 0.02 


3.8438 


0.03 ± 0.01 


0.03 ± 0.01 


2.2188 


0.08 ± 0.02 


0.08 ± 0.02 


3.0938 


0.19 ± 0.03 


0.18 ± 0.03 


3.9062 


0.01 ± 0.01 


0.01 ± 0.01 


2.2812 


0.07 ± 0.02 


0.07 ± 0.02 


3.1562 


0.03 ± 0.02 


0.03 ± 0.02 


3.9688 


0.01 ± 0.01 


0.01 ± 0.01 


2.3438 


0.08 ± 0.02 


0.07 ± 0.02 


3.2188 


0.05 ± 0.01 


0.05 ± 0.01 


4.0312 


0.03 ± 0.01 


0.03 ± 0.01 


2.4062 


0.07 ± 0.02 


0.07 ± 0.02 


3.2812 


0.03 ± 0.01 


0.03 ± 0.01 


4.0938 


0.01 ± 0.01 


0.01 ± 0.01 


2.4688 


0.08 ± 0.03 


0.08 ± 0.03 


3.3438 


0.04 ± 0.02 


0.04 ± 0.02 


4.1562 


0.02 ± 0.01 


0.02 ± 0.01 


2.5312 


0.05 ± 0.02 


0.05 ± 0.02 


3.4062 


0.05 ± 0.02 


0.05 ± 0.02 


4.2188 


0.04 ± 0.01 


0.04 ± 0.01 


2.5938 


0.07 ± 0.02 


0.07 ± 0.02 


3.4688 


0.03 ± 0.01 


0.03 ± 0.01 


4.2812 


0.00 ± 0.01 


0.00 ± 0.01 


2.6562 


0.10 ± 0.03 


0.10 ± 0.02 


3.5312 


0.05 ± 0.01 


0.05 ± 0.01 


4.3438 


0.01 ± 0.01 


0.00 ± 0.01 


2.7188 


0.13 ± 0.03 


0.12 ± 0.03 


3.5938 


0.02 ± 0.01 


0.02 ± 0.01 


4.4062 


0.00 ± 0.01 


0.00 ± 0.01 


2.7812 


0.08 ± 0.02 


0.07 ± 0.02 


3.6562 


0.02 ± 0.01 


0.02 ± 0.01 


4.4688 


0.00 ± 0.01 


0.00 ± 0.01 


2.8438 


0.10 ± 0.02 


0.09 ± 0.02 















